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ABSTRACT 

We investigate the evolution of cooling helium atmosphere white dwarfs using a full 
evolutionary code, specifically developed for following the effects of element diffusion 
and gravitational settling on white dwarf cooling. The major difference between this 
work and previous work is that we use more recent opacity data from the OPAL 
project. Since, in general, these opacities are higher than those available ten years 
ago, at a given effective temperature, convection zones go deeper than in models with 
older opacity data. Thus convective dredge-up of observationally detectable carbon in 
helium atmosphere white dwarfs can occur for thicker helium layers than found by 
Pelletier et al (1986). We find that the range of observed C to He ratios in different 
DQ white dwarfs of similar effective temperature is well explained by a range of initial 
helium layer mass between 10~ 3 and 1O _2 M0, in good agreement with stellar evolution 
theory, assuming a typical white dwarf mass of 0.6M Q . We also predict that oxygen will 
be present in DQ white dwarf atmospheres in detectable amounts if the helium layer 
mass is near the lower limit compatible with stellar evolution theory. Determination 
of the oxygen abundance has the potential of providing information on the profile of 
oxygen in the core and hence on the important 12 C(a,7) 16 reaction rate. 
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INTRODUCTION 



DQ white dwarfs have He-rich atmospheres, with the pres- 
ence of a trace of carbon. The results of model atmosphere 
calculations for 24 DQ white dwarfs have been compiled 
by Weidemann & Koester (1995). Effective temperatures, 
T eff , range from 12,500 K to 6,700 K. Log (n(C)/n(He)) 
values determined from visual spectra range from —6.2 to 
— 1.5. Analysis of IUE spectra have pushed the lower limit to 
—7.3. The generally accepted explanation for the presence 
of trace carbon in these cool He white dwarf atmospheres is 
that carbon from the core diffuses outwards to where con- 
vection dredges it to the surface. The first detailed calcula- 
tions of this process were presented by Pelletier et al. (1986) 
(hereafter P86), who show how, in the hotter phases, car- 
bon diffuses upwards from the core to be met eventually by 
the base of the surface convection zone that develops due to 
the increase in opacity following recombination of helium as 
the white dwarf cools. With further cooling more carbon dif- 
fuses upwards and the base of convection zone moves deeper 
into the star, further enriching the outer layers with carbon. 
A key finding of this research is that the maximum depth 
of the base of the convection zone nearly coincides with a 
change in ionization state of the carbon below the convec- 
tion zone. Since the slope of the carbon profile depends on 



the relative ionization states of carbon and helium, this re- 
duction in ionization state of carbon results in a steepening 
profile and a reduction in the carbon abundance in the con- 
vection zone. It is this change in ionization structure that 
causes a reduction in observed n(C)/n(He) values at lower 
temperatures, and the non-detection of carbon below 6,000 
K. Since the thickness of the helium layer regulates the effec- 
tive temperature at which the base of the convection zone 
reaches regions of diffusively enhanced carbon, the obser- 
vations of DQ white dwarfs place limits on the amount of 
helium in the white dwarf, that can be compared to the 
predictions of stellar evolution theory. P86 concluded that 
models for a white dwarf mass of M* = O.6M0 and a he- 
lium envelope mass, M env , in the range 10 -4 to 10 -3 ' 5 M» 
gave the best fit to the contemporary observational data. 
This is in conflict with evolutionary models for post AGB 
stars that predict that M env should be in the range 10 -3 
to 1O _2 M0 for this core mass. Lo complicate the picture 
further, models of the non-radial g-mode pulsations of the 
DBV star GD 358 fit the observations best for M, = O.61M 
and Af cnv = 10~ 5 7±03 A'L (Winget et al. 1994, Bradley & 
Winget 1994), which is significantly less than both the es- 
timates from the DQ white dwarf models of P86 and post 
AGB evolution models. 

In this paper, we have reinvestigated the evolution of 
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cooling helium atmosphere white dwarfs for a range of M» 
and Menv, using a full evolutionary code, specifically de- 
veloped for following the effects of element diffusion and 
gravitational settling on white dwarf cooling. The major dif- 
ference between this work and that of P86 is that here we 
use the more recent opacity data from the OPAL project 
(Iglesias & Rogers 1993). Since, in general, these opacities 
are higher than those available ten years ago, at a given 
effective temperature, convection zones go deeper than in 
models with older opacity data. Thus convective dredge-up 
of observationally detectable carbon in helium atmosphere 
white dwarfs can occur for thicker helium layers than found 
by P86. We find that the range of observed C to He ratios 
in different DQ white dwarfs of similar effective tempera- 
ture is well explained by a range of M cnv between 1CF 3 and 
1O _2 M0, in good agreement with stellar evolution theory, 
assuming a typical white dwarf mass of 0.6M Q . 

The details of the evolutionary code are given in the 
next section. In section 3, we present our results for the 
evolution of the abundance profiles and compare them with 
observations and current ideas for the evolution of possible 
white dwarf progenitors. Our conclusions are given in section 
4. 

2 MODELING DIFFUSION AND 
CONVECTIVE MIXING IN COOLING WHITE 
DWARFS 

To study the evolution of the carbon distribution in helium 
atmosphere white dwarfs, we have added routines that solve 
the equations describing element diffusion and gravitational 
settling to a stellar evolution code specifically designed for 
following the cooling history of white dwarfs. We restrict 
attention to mixtures of helium, carbon and oxygen. The 
existence of OPAL opacities for arbitrary mixtures of these 
three elements allows us to calculate consistently the effects 
of element diffusion on the optical thickness of the non- 
degenerate envelope, which controls the rate of cooling of 
the white dwarf. These effects were not included by P86. 
The initial model consists of a core of uniform composi- 
tion, surrounded by an envelope also of uniform composi- 
tion. The initial temperature profile is determined by com- 
puting a thermal equilibrium model, in which radiative and 
neutrino losses are balanced by an artificial constant energy 
rate, e c . We study three sets of models, labeled HE, LO, 
and PG. To identify our sequences, we use an alphanumeric 
label that consists of three parts. The first part identifies 
the core mass, the second part the particular set of models 
to which the sequence belongs, and the third part identi- 
fies the envelope mass. For example, the P6HE6 sequence is 
part of the HE set of models and has a core mass of 0.6Mq 
and envelope mass of approximately 1O _6 M0. The HE set 
is the largest set of models and explores how the degree of 
carbon and oxygen dredge-up depends on core mass and en- 
velope mass. The LO set has a lower core oxygen abundance 
than the HE set and is used to explore the relation between 
the oxygen abundance in DQ photospheres and the oxygen 
abundance in the core. The PG sequence has envelope com- 
position characteristic of the PG1159 stars and is used to 
test the idea that some DQ white dwarfs are descendants of 
the PG1159 stars. In particular, Dehner & Kawaler (1995) 
have suggested that the unusual DBQv GD 358 is such an 



object. The properties of our initial models are given in Ta- 
ble 1. M corc , Xc.corc, Ao.coro and M env , Xc,cnv, ^o.cnv are 
the initial masses, carbon mass fraction and oxygen mass 
fraction of the CO core and He-rich envelope respectively. 
Our range of M cnv values at O.GMq was chosen to include 
those from the asteroseismological studies of GD 358 up to 
the largest values consistent with stellar evolution models 
for the white dwarf progenitor. 

We stop our evolutionary calculations when the temper- 
ature at zero optical depth is 6,000 K, which is the lowest 
temperature in the OPAL opacity tables. This corresponds 
to T cB = 7, 000 K. 

The treatment of diffusion is fairly straightforward. We 
use an approach that is applicable to a general mixture of an 
arbitrary number of elements. We neglect thermal diffusion, 
which has been shown to be negligible for white dwarfs by 
Iben & MacDonald (1985), and Paquette et al. (1986b) and 
radiative levitation. Although radiative levitation is impor- 
tant for determining the photospheric composition of hot 
white dwarfs (Fontaine & Michaud 1979, Vauclair, Vauclair 
& Greenstein 1979), it is entirely negligible in the interior 
and also in the photospheric regions at the effective temper- 
atures of the DQ white dwarfs. For simplicity, we also ne- 
glect magnetic fields and stellar rotation. The composition 
profiles are then determined by the competition between 
gravity, partial pressure gradients and induced electric fields. 
The diffusion velocities satisfy (Curtiss & Hirschfelder 1949, 
Aller & Chapman 1960, Burgers 1969, Muchmore 1984, Iben 
& MacDonald 1985) the multi-fluid equations 

^1 _ £l± _ niZi eE = V Ka (w 3 - Wi ) (1) 
dr p dr ' 

3 

where p», pi, m, Zi, Wi are the partial pressure, mass den- 
sity, number density, mean charge, and diffusion velocity for 
species i and E is the electric field induced by the gradi- 
ents in the ion densities. The resistance coefficients are from 
Paquette at al (1986a). 

In solving equations (1) for the diffusion velocities, care 
must be taken to avoid rounding errors due to the subtrac- 
tions. This is particularly important for the diffusion veloc- 
ities of trace elements. For N species including electrons, 
equations (1) is a set of N — 1 independent linear equations 
for N + 1 unknown variables (the Wi's and E). To complete 
the set of equations, we use the conditions for no net mass 
flow relative to the center of mass 

^2 l A i n i w i =Q (2) 

i 

and charge neutrality 

Z i n i = ( 3 ) 

i 

In equations (2) and (3), the sums are over all species (ions, 
neutrals and electrons). We assume the electron to have zero 
mass. The electrical field can then be eliminated directly 
since equation (1) applied to electrons gives 
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The electron pressure, p e , and the electron number density, 
n e , can now be eliminated since 



Pe =p-^Pi 
and 



= ^2 Z ^ 



(5) 



(6) 



where the sums now do not include the electrons. Using 
the equation of hydrostatic balance, equation (1) can be 
transformed to 



dpi _ x iijZj 
dp 1 Y. n i Z 3 



E 



dp 



(7) 



To reduce rounding errors, we reformulate equation (7) as 



E A 

dpi _ x \ 3^i,3^e 1 

dp J X A Zj 



E O .v) 



E x i z i 

, A 3 
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The equations for the changes in mass fractions due to ele- 
ment diffusion and convective mixing are 



dXi 
~dT 



OF, 







dm 

where the fluxes, Fi, are 



Fi = 4nr 2 pWiXi + a con ^l 
om 



(9) 



(10) 



Convective mixing is treated as a diffusion process, with 
diffusion coefficient, a con = (Anr 2 pi) 2 /t CO h where I is the 
mixing length and r CO n is the convective turnover time scale 
calculated from mixing length theory. 

Equations (9) and (10) for each element are solved to- 
gether with the equations of stellar structure by standard 
implicit finite difference techniques. The diffusion velocities, 
Wi, are first found by standard matrix methods (LU de- 
composition), from finite difference equations derived from 
equations (8). 

This multi-fluid formulation requires knowledge of the 
partial pressures, pi. For non-interacting species or species 
interacting through short-range forces such as Van der Waals 
forces, calculation of pi is relatively straightforward. How- 
ever for species interacting through long-range coulomb 
forces, it is not clear that the concept of partial pressures 
is meaningful. We circumvent this difficulty by adopting a 
'linear mixing law', 



E: 



' E n i 



(11) 



3 EVOLUTION OF THE ABUNDANCE 
PROFILES 

Figure 1 shows how the carbon mass fraction changes with 
T e ff for the P6HE3 sequence. Also shown are the location of 
the photosphere and the boundaries of the convection zones. 
This figure can be compared to figure 1 of P86. The most im- 
portant difference is in M cmax , the maximum depth of the in- 
ner boundary of the convection zone. For M» = 0.6M Q , and 
M onv = 1.9 10~ 4 M Q , P86 find M cmax = 3.8 10" 7 M Q , which 
occurs when T off is near 10 4 K. For the P6HE4 and P6HE3 
sequences, we find M cmax = 1.8 10~ 6 and 1.7 10~ 5 M Q , 
respectively. We attribute most of the order of magnitude 
difference in M cmax to differences in radiative opacity. Al- 
though the depth of the convection zone is sensitive to con- 
vective efficiency for relatively high effective temperatures, 
Fontaine, Tassoul & Wesemael (1984) have shown that at 
temperatures characteristic of DQs the dependence is sig- 
nificantly less. To illustrate this, we have calculated cool- 
ing sequences without element diffusion for a 0.61Mq white 
dwarf with M anv = 10 _ M@ and mixing length ratios, a — 
1.0, 1.5 and 2.0. In figure 2, we show the location of the 
convection zone inner boundary for pure helium envelopes 
and envelopes that are a mixture of helium and carbon with 
Ac = 10~ 3 . To test the dependence of M cm ax on opacity, we 
have also calculated cooling sequences using the Los Alamos 
opacities (Huebner et al. 1977). In figure 3, we show the lo- 
cation of the convection zone inner boundary for the helium 
and carbon mixture for the OPAL and Los Alamos opac- 
ities. It is clear that the newer opacities result in deeper 
convection zones and that the depth of the convection zone 
is insensitive to a for T c g < 14, 000 K. We also see that the 
presence of carbon in small amounts reduces the depth of 
the convection zone, and hence the convection zone in the 
P6HE4 sequence is smaller than that of the P6HE3 sequence. 

3.1 Comparison with observations and prior 
evolution 

In figure 4, we show the ratio of the number densities of 
carbon and helium in the convection zone as a function of 
T e ff for the models with M CO rc = 0.6M©. Also shown are 
the observed photospheric ratios taken from Weidemann & 
Koester (1995). When possible, error bars are from the in- 
dividual papers referenced in Weidemann & Koester (1995). 
Excluding the stars for which there are only upper limits, 
the mean of the observed C to He ratios is log n(C)/n(He) 
= —5.35 and the corresponding mean T e s is 9,300 K. For 
our models with M cor o = 0.6M Q , these values are obtained 
if M onv = 1.1 1O" 2 M . 

There are many possible ways that helium atmosphere 
white dwarfs can form. They may be descendants of asymp- 
totic giant branch (AGB) stars or sdOs. They may also be 
the result of mergers in double degenerate systems. Because 
the evolutionary state of sdO stars is poorly understood and 
the physics of mergers is very complex, we assume here that 
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Table 1. Properties of the initial models 



Sequence 




^C,corc 


^O,core 


M env 




^0,env 


R 




Pc 


T c 




U.OUU 


U.OU 


U.OU 


4.yo X 1U 


U.Ul 


U.Ul 




1 .40 


z. f o 


1 f!7 

l.U/ 


P6HE4 


0.600 


0.50 


0.50 


9.36 x 10" s 


0.01 


0.01 


1.16 


7.40 


2.78 


1.07 


P6HE3 


0.600 


0.50 


0.50 


9.60 x 10~ 4 


0.01 


0.01 


1.14 


7.45 


2.80 


1.07 


P6HE2 


0.602 


0.50 


0.50 


8.29 x 10~ 3 


0.01 


0.01 


1.13 


7.43 


2.97 


1.06 


P6L03 


0.600 


0.99 


0.01 


1.02 x 10~ 3 


0.01 


0.01 


1.16 


7.49 


2.69 


1.03 


P6L02 


0.602 


0.99 


0.01 


8.29 x 10~ 3 


0.01 


0.01 


1.15 


7.47 


2.86 


1.07 


P6PG6 


0.600 


0.50 


0.50 


4.41 x 10~ 6 


0.50 


0.17 


1.88 


10.97 


1.40 


2.02 


1HE4 


1.00 


0.50 


0.50 


1.01 x 10~ 4 


0.01 


0.01 


0.590 


7.68 


30.7 


1.08 


1HE3 


1.00 


0.50 


0.50 


1.02 x 10~ 3 


0.01 


0.01 


0.590 


7.73 


30.9 


1.10 


Units are 


Mq for 


masses, 10 9 


cm for R, 10 4 K for T eS , 


10 6 g cm" 


~ 3 for p c 


, and 10 8 K for T c 
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-2 1 1 ' 1 ' 1 1 1 1 1 1 1 ' 1 ' 1 ' 1 ' 1 1 
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loot,, 



Figure 1. Convection zone boundaries (solid line), location of the 
photosphere (dotted line) and contours of carbon mass fraction 
(dashed lines) for the P6HE3 sequence. The contour levels are, 
from top to bottom, log Xci2 = —10, —9, —1. 
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4.8 4.7 4.6 4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 
LogT„(K) 



Figure 2. Location of the inner boundary of the convection zone 
for a 0.61Mq white dwarf with M env = 10 -2 Mq and mixing 
length ratios, a = 1.0, 1.5 and 2.0 (from bottom to top). The solid 
lines are for Xq = 10~ 3 and the dashed lines are for Xc = 0. 



DQs have AGB predecessors. 

The possible modes of how AGB stars evolve to white 
dwarfs have been described in detail by Iben (1984). In con- 
trast to the case for DA white dwarfs, there have been too 
few detailed calculations of the evolution to helium atmo- 
sphere white dwarfs for us to be precise about the mass of 
helium in these objects. However due to the relationship be- 
tween the white dwarf and its AGB predecessor, we find it 
convenient to parameterize M cnv in terms of M^!* , the mass 
of the helium layer just prior to ignition of a helium shell 
flash. This is given by (Iben & Tutukov 1996) 

log Mh c " = -1.835 + 1.73M core - 2.67M c 2 oro (12) 

for Mcorc > O.55M . This is in good agreement with the 
results for the final helium layer mass from unpublished cal- 
culations by one of us (JM) of the evolution of population 
I helium stars of masses between 0.5 and 0.8 Mq to the 
white dwarf stage. For less massive cores, we find that this 
equation underestimates Mh" 1 '. 

Whether the final white dwarf has a hydrogen-deficient 
envelope or not depends on the phase of the thermal pulse 
cycle at which the star leaves the AGB. If this occurs dur- 
ing helium shell burning, winds can remove all remaining 



hydrogen. Thus we might expect M env to be close to Mg" , 
the mass of the helium layer at the end of quiescent he- 
lium burning, in a star that is near the tip of the AGB 
(D'Antona & Mazzitelli 1991). This is about 30% of Mg* 
(Vassiliadis & Wood 1993; Iben & Tutukov 1996). D'Antona 
& Mazzitelli (1991) estimate that M cnv will be about 30% 
less than M^" d because carbon and other heavy elements 
will sink out of the helium buffer layer as the star evolves 
to a white dwarf. However, Iben & Tutukov (1984) find 
for M» = 0.6M Q , that M onv = 0.016M Q , which is about 
0.9 M^j*. Thus Menv/M^j* is of order 0.2 to 1.0 so that 
for M C ore = O.6M0, we estimate M cnv is between 3.5 10~ 3 
and 1.7 1O _2 M0, which is in reasonably agreement with the 
mean value derived above for the DQ white dwarfs. Fur- 
thermore this range of M cnv corresponds to a range in log 
n(C)/n(Hc) at 9,300 K of -4.2 to -5.8. Although many of 
the DQs lie within these limits, clearly some of the data 
points fall outside this range. The DQ white dwarf with the 
largest C to He ratio, G35-26, has log n(C)/n(He) = -1.5 
and T off w 12, 500 K. If M coro = 0.6M Q , this would require 
M onv = 2 1O~ 5 M0. At the other extreme, there are two DZ, 
two DB and 1 DBA white dwarfs with upper limits to log 
n(C)/n(He) of less than —6.5. 

We suggest that DQs with very high C to He ratios are 
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Figure 3. Location of the inner boundary of the convection zone 
for a 0.61-Mq white dwarf with M en v = 10 — 2 Mq and mixing 
length ratios, a = 1.0, 1.5 and 2.0 (from bottom to top). The 
solid (dashed) lines are for when OPAL (Los Alamos) opacities 
arc used. 



freB charges. We use the occupancy probability formalism 
of iummer and Mihalas (1988). Fontaine et al (1977) di- 
vi( i the temperature- density plane into three regions. In 
a w density regime, they truncate the partition function 
by iumming over the finite number of bound state energy 
le-\ Is calculated from the static screened coulomb poten- 
tia In a high density regime they assume complete ioniza- 
tic . In the intermediate region they interpolate between 
th low and high density regimes. In table 2, we give some 
rel vant properties of the convection zones for the coolest 
m( lei of each of the HE sequences. In each case, conditions 
an such that they lie in the high density fully ionized regime 
of he Fontaine et al (1977) equation of state. We have ai- 
re; ly noted that P86 find shallower convection zones. At a 
s depth of 4.0 1O~ 7 M0, corresponding to the depth of 
base of the convection zone in the model of P86 with 
= 1.9 1CT 4 M Q , the temperature is 1.7 10 6 K and the 



sity is 7.7 10 1 g cm 3 for the P6HE4 model, which again 



m< 
th 
M 
de 

places this point in the high temperature fully ionized region 
of the Fontaine et al (1977) equation of state. Hence we con- 
clude that the difference in ionization structure in the layers 
neighboring the base of the convection zone must be due to 
higher temperatures in our models than those of P86. 



significantly more massive than the typical white dwarf. Due 
to the higher gravity these have smaller scale heights in the 
outer layers, which leads to less massive helium envelopes. 
Figure 5 is the same as figure 4 except that the core mass 
is 1.0 Mq. The solid lines are for envelope masses of 10~ 4 
and 10 _3 Mq. The range in M cnv consistent with prior evo- 
lution is 5.0 10~ 4 < Af onv < 1.7 1O~ 3 M . Extrapolating our 
results at 0.6 and 1.0 Mq, we estimate that the minimum 
core mass consistent with the observed C to He ratio for 
G35-26 is about 1.3 Mq, which is consistent with the high 
gravity found by Thejll et al. (1990) for this star from anal- 
ysis of spectra. White dwarfs of mass greater than 1.2 Mq, 
however, are thought to evolve from stars of initial mass 
between 7.5 and 10 Mq (with an uncertainty of about 0.5 
Mq in both limits) which experience quiescent core carbon 
burning and so develop ONe cores (Dommguez, Tornambe 
& Isern, 1993; Ritossa, Garcfa-Berro & Iben 1996). Hence 
we can only conclude that the white dwarf in G35-26 is more 
massive than 1 Mq. 

Similarly the lowest C to He ratios can be explained 
by allowing the white dwarf to be less massive than the 
canonical 0.6 Mq. The maximum M cnv for a 0.5 Mq white 
dwarf is 0.04 Mq. With this much helium, we estimate that 
n(C)/n(He)~ 10~ 7 , which is about the lowest detected car- 
bon abundance. 

In contrast to P86, we do not find that the carbon abun- 
dance starts to drop near T e s — 7, 000 K. This difference 
is due to carbon being fully ionized below the convection 
zone in the P6HE4, P6HE3 and P6HE2 sequences, whereas 
P86 find partial recombination. In their calculations, P86 
used the Fontaine, Graboske and Van Horn (1977) equa- 
tion of state. Our equation of state is similar in that it is 
also based on free energy minimization and is thermody- 
namically consistent. An important difference that might be 
the cause of the difference in ionization structure is in the 
treatment of Stark ionization i.e. the dissolution of bound 
states due to the ambient electric field from the presence of 



3.2 The role of oxygen 

The mass fraction of oxygen in the white dwarf core has a 
strong dependence on the rate of the 12 C(q, 7) 16 reaction. 
Although recent experimental progress has reduced the error 
in this reaction rate to about 35%, this uncertainty can still 
be significant for nucleosynthesis studies and models of type 
I supernovae (Arnett 1996). Hence a measurement of the 
oxygen abundance in the atmosphere of a DQ white dwarf 
may provide useful information about this reaction rate. 

Salaris et al. (1997) have recently calculated the evo- 
lution of population I stars of mass 3.2 and 7.0 Mq, using 
the best current estimate for the 12 C(a, 7) 16 reaction rate 
(Woosley, Timmes & Weaver 1993; Thielemann, Nomoto 
& Hashimoto 1996; Arnett 1996). These stars develop CO 
cores of mass 0.6 and 1.0 Mq. The central mass fractions are 
Ac = 0.233, X = 0.739 and X c = 0.316, X = 0.655 re- 
spectively. In our unpublished calculations of the evolution 
of population I helium stars, using the Arnett (1996) reac- 
tion rate, we find slightly higher central oxygen abundances. 
In the context of dredge- up of oxygen diffusing from the core 
a more relevant quantity is the O mass fraction just below 
the extinct helium burning shell. The Salaris et al. (1997) 
models give Xo w 0.02 at the center of the helium shell. 
For the population I helium stars, we find Xo ~ 0.05 at the 
same location. 

To get some idea of how the abundance of oxygen 
in the core affects the amount of oxygen dredged up to 
the surface, we have made two additional evolutionary se- 
quences, P6L03 and P6L02 for 0.6 Mq cores of composi- 
tion Xc = 0.99, Xo — 0.01. Although the abundance pro- 
files in our white dwarf models differ in detail from those 
of the evolutionary calculations, the two values of core oxy- 
gen abundance considered here roughly bracket the possible 
range of oxygen abundance in the region below the helium 
shell. In figure 6, we show the O to He ratio as a function 
of T e ff for these two sequences and the P6HE2 and P6HE3 
sequences. In the sequences with the thicker helium layer 
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Table 2. Properties of the convective zones for the 
coolest model of each HE sequence 



Sequence 








'Tconv 


pconv 


P6HE6 


9,930 


1.973 x 10" 


-7 


3.08 


33.7 


P6HE4 


7, 640 


6.950 x 10- 


-7 


2.26 


117 


P6HE3 


6,830 


7.015 x 10" 


-6 


1.07 


740 


P6HE2 


6,960 


7.621 x 10" 


-5 


0.857 


3610 



Units are Mq for masses, K for T eff , 10 6 K for T c 
and 10 2 g Cm for Pconv • 




Figure 4. The solid lines are the values of log ra(C)/n(He) for 
the P6HE6, P6HE4, P6HE3 and P6HE2 sequences. Also shown 
with error bars arc the observed ratios. 




Figure 5. As figure 4 but for the 1HE4 and 1HE3 sequences. 

very little oxygen is dredged to the surface, independent of 
the amount of oxygen in the core. For both of the thinner 
helium envelopes, oxygen is expected to be present in de- 
tectable amounts (based on a comparison of the oscillator 
strengths of the O I multiplets at 130.5 and 777.4 nm with 
the detected C I lines). Also the maximum photospheric 
oxygen abundance depends on the core oxygen abundance. 
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he P6HE3 sequence with Xo.core = 0.5, the maximum 
tospheric abundance is log n(0)/n(He) = —4.0 which 
irs at T efI = 8.0 10 3 K and in the P6L03 sequence with 
o = 0.01, the maximum photospheric abundance is 
n(0)/n(He) = -5.3 which occurs at T cB = 8.3 10 3 K. 



He ce measurement of the photospheric oxygen abundance 
in , DQ has the potential of providing useful information 
ab lit the interior composition of the white dwarf. 



The DBV GD 358 

U 1 spectroscopy of the DBV GD 358 by Provencal et al. 
(l! )6) has confirmed the detection of He II 164.0 nm and C 
II 
ca 
27 
al 
he 



th 
an 
He 

to 
st« 
ag 
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53.5 nm, first reported by Sion et al. (1988). Thus GD 358 
be considered to be the hottest DQ star, with T e ff « 
)00 K and log n(C)/n(He) = -5.65 ±0.10 (Provencal et 
996) . At this temperature, the star is sufficiently old that 
ry elements will have sunk deep beneath the photosphere 
yet too hot for a deep convection zone to develop and dredge- 
up carbon diffusing out of the core, unless the helium layer 
is very thin. 

Dehner & Kawaler (1995) have investigated the possibil- 
ity^hat GD 358 is a descendant of PG1159. They calculate 
evolution of the abundance profiles in a star that has 
nitial envelope of mass 3 1O _3 M0 and composition 30% 
35% C and 35% O by mass. When the star has cooled 
he temperature range of the DBV's, they find that the 
has a layer of pure helium of mass 1O _5 ' 5 M0, in good 
ement with the asteroseismological analysis of the WET 
a-vations by Bradley & Winget (1994). Our P6PG6 se- 
qu ice begins with a helium mass comparable to that found 
by isteroseismology for GD 358. When T cS = 27, 000 K, 
th center of the transition from He to C is at a mass depth 
of 0~ 6,7 Mq, i.e. an order of magnitude less than found 
frc l asteroseismology, yet the photospheric n(C)/n(He) is 
0" 20 , which any orders of magnitude smaller than 
ob ?rved. To get a detectable C to He ratio at 27,000 K, 
th helium layer mass cannot be too much larger than the 
section zone, which at this temperature has its base at 
m; s depth of 1O~ 14 M0. Hence scenarios in which GD 358 
is descendant of the PG1159 star suffer from a conflict 
between the requirement that the helium layer mass be con- 
sistent with asteroseismology and the requirement that it 
be very small for C to be present in the photosphere at the 
observed abundance. 

Other possible sources for the photospheric carbon are: 
1) that it is radiatively supported, 2) accretion from the 
ISM and 3) that it is left over from a merger of two low 
mass white dwarfs in a double degenerate binary. 

Detailed calculations of radiative levitation in helium 
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Figure 6. The values of log n(0)/n(He) for the P6HE3 and 
P6HE2 (solid lines) and the P6L03 and P6L02 (dashed lines) 
sequences. 

dominated atmospheres have recently been presented by 
Chayer, Fontaine & Wesemael (1995), who find that, for 
T c s < 4 10 4 K, the radiatively levitated abundance of carbon 
is significantly less than n(C)/n(He) = 10~ 6 . Hence radia- 
tive levitation appears to be insufficient to support carbon 
at the levels observed in GD 358. 

Clearly, if accretion is invoked as the mechanism for the 
presence of photospheric carbon, then an explanation of why 
hydrogen has a low photospheric abundance of n(H)/n(He) 
< 3 10 (Provencal et al. 1996) is needed. Mechanisms 
that prevent or inhibit accretion have been reviewed by Al- 
cock & Illarionov (1980). However most of these mechanisms 
work equally well for carbon as for hydrogen. One mecha- 
nism that permits differential accretion is that the white 
dwarf has sufficiently high rotation rate and magnetic field 
that the 'propeller' mechanism (Illarionov & Sunyaev 1975) 
reduces the accretion rate of ionized hydrogen below the 
Bondi-Hoyle rate. A significant fraction of the carbon in 
the interstellar medium (ISM) can be in the form of grains, 
which have a significantly lower charge to mass ratio than 
protons so that they follow independent particle trajecto- 
ries while outside the white dwarf's magnetosphere (Alcock 
& Illarionov 1980). If they cross into the magnetosphere and 
are then evaporated by the radiation from the white dwarf, 
they provide a source of ionized carbon and other species 
that can then follow magnetic field lines to the white dwarf's 
surface. The requirements for an effective propeller are dis- 
cussed by Alcock & Illarionov (1980) and Stella, White & 
Rosner (1986). Accretion of ionized hydrogen is centrifu- 
gally inhibited if the co-rotation radius, r c , is less than the 
magnetospheric radius, r m . By studying multiplet splitting 
in their WET observations of the oscillations of GD 358, 
Winget et al. (1994) find evidence for a magnetic field of 
average strength 1300± 300 G and for differential rotation 
with the outer envelope rotating with a period of 0.89 d. 
With these parameters, r c w 5 10 9 cm. To determine r m , 
we equate the magnetic energy density to the kinetic energy 
density on the symmetry axis of the (upstream) accretion 
flow. Assuming a dipole magnetic field of strength _B*, 
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wl re Voo and rice are the velocity of the ISM relative to 
st< and the ISM number density respectively. For GD 358, 
Pr vencal et al. (1996) find « 40 km s _1 and hence 
> r c if iioo < 3 10 5 cm~ 3 , which is the case for much of 
volume of the ISM. Hence the propeller mechanism is ca- 
le of inhibiting accretion of ionized material by GD 358 
this may be the reason why the photospheric abundance 
ydrogen is low. However, for carbon to be accreted at a 
gn iter rate than hydrogen, we also require that the charac- 
tei itic distance at which incoming grains evaporate, r c , be 
than r m - The temperature of dust grains of character- 
size 0.3 fim heated by the radiative flux from a star of 



distance r and luminosity L is (Osterbrock 1989) 
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For refractory grains which evaporate at « 10 K, r c w 10 ' 5 
cm which is much greater than r m for any values of n rel- 
evant to the ISM. Hence it is unlikely that the propeller 
mechanism can inhibit accretion of hydrogen without doing 
the same for carbon. From our P6PG6 sequence, we estimate 
that the time scale for settling of carbon out of the convec- 
tion zone is about 0.1 yr. Hence, in a steady state, carbon 
has to be accreted at a rate of 4 1O _19 M0 yr -1 to maintain 
the observed photospheric abundance. For fluid dynamical 
accretion, and cosmic abundances this requires an ISM den- 
sity of riH = 5 cm -3 , which is higher than in the local ISM. 

The merging of two white dwarfs in a double degener- 
ate scenario is a complex process that is poorly understood. 
However it might reasonably be expected that the result of 
the merger is a luminous object that has much evolution 
to go through before reaching properties similar to those of 
GD 358. Hence given the short time scale for the evolution 
of the carbon distribution, it would be a remarkable coin- 
cidence that we find the star just at the time it looks like 
GD 358. 

Hence none of these mechanisms can satisfactorily ex- 
plain the observed photospheric abundances in a way that 
is consistent with the data from asteroseismology. 



4 CONCLUSIONS AND DISCUSSION 

Our main conclusion is that use of up-to-date opacities in 
modeling the evolution of DQ white dwarfs removes the dis- 
crepancy between estimates of the helium layer masses from 
the observed carbon abundances and the predictions of stel- 
lar evolution theory. In particular the bulk of the DQ white 
dwarfs can be interpreted as having mass O.6M0 and helium 
layer masses between 10 -3 and 10 _2 Mq. 

We suggest that DQ white dwarfs with exceptionally 
high atmospheric C abundances, such as G35-26, are of 
higher than average mass and those with no C or very low 
C have lower than average mass. 
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We also predict that oxygen will be present in DQ white 
dwarf atmospheres in detectable amounts if the helium layer 
mass is near the lower limit compatible with stellar evolu- 
tion theory. Determination of the oxygen abundance has the 
potential of providing information on the profile of oxygen in 
the core and hence on the important 12 C(a,7) le O reaction 
rate. 

In contrast to P86, we do not find a drop in photospheric 
carbon abundance near T e g — 7,000 K. Unfortunately, we 
are not able to continue our calculations below this tempera- 
ture due to lack of opacity tables for low temperatures. Fur- 
thermore, for T e ff below about 10,000 K, uncertainties in 
the equation of state are becoming increasingly important. 
In particular, the equation of state is becoming sensitive to 
the treatment of Stark ionization. Hence we do not rule out 
the possibility that partial recombination of carbon occurs 
below the convection zone leading to a drop in photospheric 
carbon abundance as found by P86. 
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